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CHAPTER 1 : INTRODUCTION 


Determining the crystal orientation by the Laue 
method involves indexing some of the Spots and/or some of the 
Zones of the Spots on the Laue photograph. Indexing a crystal 
which is oriented at random with respect to the X-Ray beam 
proves difficult. The routine for indexing Laue back 
reflection patterns can be tedious and time consuming. Even 
for a skilled operator, conversant with Greninger chart, 
Wulff net and Stereographic projections, it often requires 
considerable time and experience before difficult patterns 
can be solved. Indexing of Laue patterns for cubic crystals 
is generally easier, since the tables of interplanar angles 
and standard stereographic projections are available. When 
the crystal structure is non-cubic, tables of interplanar 
angles may not be available and indexing a pattern is then 
not possible. Since the labour involved in indexing Laue 
patterns is enormous, various approches have been made to 
index the pattern by computer - 

The present work involves developing a 
computer algorithm for the rapid solution of Laue back 
reflection patterns of body centered tetragonal crystals. 
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Reflections of planes having indices (666) to (66 6 ) 

are indexed directly by using the X and Y coordinates of 
the r>pots on the film as the input data. The code developed 

finds out the indices of the important zone axes and the 

mutual angles of these axes are also matched with the 
given angles obtained from the stereographic projection of 

the pattern. The computational scheme compares the 

interplanar angles for the various spots on the film with the 
theoretical interplanar angles, which are calculated by the 
program after assigning Miller indices to the planes causing 
these spots. The systematic error associated with measurement 
of the distance of separation between the seimple and the 
X-Ray film is minimised to determine the best value of 
distance. 

Directionally solidified alloys of Tin have 

dendritic structure. The primary dendrites are plate like in 
shape and are roughly parallel to the growth direction, which 
in turn is the ingot axis. It is noted from the literature 
that, the growth direction in these alloys is [ 1 1 0 3 

This value has been reported for pure Tin [1], as well as for 
Sn-12% Bi alloys [2,33- The present work involves developing 
the computer technique for indexing the LBR patterns of these 
alloys and also to check the obtained crystallographic 
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orientation with the pubJished data. 

Tin alloys with 2 %, 6% and 12% Bismuth were used for 
obtaining the Laue back reflection patterns for subsequent 
indexing. These alloys are directionally solidified at 
different freezing rates. The crystallographic orientation of 
each sample is determined with the help the indices of the 
zone axes obtained from the program. This is then compared 
with the published data for the orientation of Sn-Bi alloys. 

The computer program entitled " LAUE " , is written 
in FORTRAN 77 and the executions are carried out in Hewelett 
Packard 9000 system. The execution takes about 2 to 3 seconds 
of CPU time for indexing 7 spots, with an allowable error of 
3.0 degrees. 
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CHAPTER 2 : LITERATURE SURVEY 


The various methods of indexing Laue patterns using 
computer are briefly reviewed in this chapter. 

A number of workers have tried indexing the Laue 
patterns by computer. Many of the techniques adopted are for 
a specific crystal system and the others are for any crystal 
system. One of the early attempts is by D.T.Camp and J.A.Clum 
[4]. They have developed a scheme for indexing LBR patterns 
of Zinc crystal (hexagonal). The scheme first calculates the 
interplanar angles between the planes with indices, 
within a specified limit (333) to (333 ) and stores 
them. These angles are then compared with the experimental 
interplanar angles measured from a stereographic projection. 
The measured input data is systematically compared with the 
generated table of steindard angles. All combinations of poles 
whose interplanar angles compare with the requisite specified 
accuracy, are printed out. The computer program was written 
in FORTRAN IV. 


Christensen et al. [5,6] have developed a technique 
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which dupJi cates the traditional technique of matching the 
angles between prominent reflections on the film, with the 
true values of interplanar angles. The scheme seeks a 
solution by assigning trial indices to a pair of reflections 
generating the theoretical interplanar angle and comparing it 
with its counterpart from the X-Ray pattern. They have 
limi.ed their search between (333) to (33 3 ), for 
the planes causing the spots. They have adopted this scheme 
for indexing Silver single crystals (cubic). The same 
technique has been used by Cyril Anasia et al. [7] for 
indexing Sapphire single crystals (hexagonal). The 
computer program was written in FORTRAN IV. 

Haskell V. Hart et al. [8] have developed a scheme 
of indexing LBR patterns of echinoid Calcite (cubic). Here 
they identify the spots belonging to each zone and they 
assign indices to these zone axes. The angles between zone 
axes are measured and compared with the experimentally 
calculated angles, (using the r and & values of the 
hyperbolas from the Greninger’s chart). They have written the 
computer program in BASIC. A similar approch of indexing the 
zone axes has been adopted by Fewster [9] also. 


Jean Laugier et al. [10] have developed a 


scheme 
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which does •the indexing of LBR patterns of any crystal 
system. The computer program was written in FORTRAN IV. The 
data observed ar<i ( X,Y ) co-ordinates of the spots on the 
Laue photograph. Two reflections which are at the 
intersection of zonal planes and thus, having small Miller 
indices are indexed first. The angle between the 
corresponding scattering vectors is compared, within a given 
angular tolerance, to the angles between all possible pairs 
of reciprocai vectors upto a given limit for the Miller 
indices. The verification of calculated solutions is made by 
means of the visual comparison between the experimental Laue 
photograph and the simulated patterns displayed on the same 
scale. Thus, by simple super imposition of the observed and 
the calculated patterns, all the observed spots are checked 
visually. 


The computer program developed by C.A. Cornelius 
[11] also involves comparison of the Laue pattern generated 
by computer with the experimentally obtained photographs. 
This scheme is for indexing LBR patterns of any crystal 
structure . 

The present work incorporates the scheme similar to 
the one developed by Christensen et al. [5,6], which requires 
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measuring the X and Y co-ordinates of the spots on the film. 
In addition to matching the interplanar angles, the zone axes 
are also indexed and the mutual angles between these axes are 
matched with the values obtained from the stereographic 
projection of the pattern. 

The present work also involves determining the 
crystallographic orientation of Sn-Bi alloys and comparing it 
with published data. The primary dendrites of directionally 
solidified alloys of Sn are plate like in shape and are 
roughly parallel to the growth direction. G.F. Bolling et al . 
[1} have determined the growth direction of pure Tin (99.989 
+ purity) to be I 110 ]. P.J. Aliearn et al. [2] have 
confirmed this growth direction for Sn-12% Bi alloy also. 
J.C. Warner et al. [3] have also reported f 1 1 0 ] as the 
growth direction for Tin rich alloys. 

The computational scheme involves calculation of 
theoretical interplanar and interzonal angles. These values 
have been checked with published values for Tin crystals by 
R.E. Trounfelker et al. [13] and B.S. Chandrasekar et al. 
[14]. The lattice constants c and a are obtained from the 
work published by J.A. Lee et al, [15]. Lee et al. report a 
change in c/a ratio from 0.5455 for pure tin to 0.5451 for 6% 
Bi in Tin. The effect of such a variation in C/a ratio on an 
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interplanar angle of 28.8298 degree.s (between (12 3) and 
(332)) is found to be 0.007A2 degrees. Since this 
variation in interplanar angle due to a change in c/a ratio 
is negligible compared to the permissible error of 3.0 
degrees used in the computer code, it was decided to use C/a 
ratio corresponding to pure Tin itself, which is 0.5455. 
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CtlAPTER 3; EXPRIMENTAL WORK 


The experimental work carried out in this project 

involves, 

(i) preperation of the sample for X-Ray exposure, 

(ii) obtaining the Laue back reflection pattern from 
the sample, 

(iii) obtaining the input data from these LBR patterns. 


3 . 1 Preperation of the samples : 

Directionally solidified Tin - Bismuth samples were 
used for taking the Laue back reflection patterns. 
Directionally solidified ingots of Tin with 2%, 6% and 12% 
Bismuth had been prepared by Mr. N.Das in the Solidification 
Laborotory earlier. The ingots were solidified at different 
freezing rates. TABLE I shows the freezing rate for different 


compositions . 
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TABLE 1. EXPERIMENTAL SAMPLES 


Wt % Bi 

Rate of 
freezing U'/s) 

Sample No. 

Computer output 
No. (Appendix III) 

2 

41.7 

4.3 

5 


68.4 

2.1 

3 


255.0 

3.1 

4 

6 

41.7 

6.4 

10 


41.7 

6.4 M 

8 


68.4 

6.2 

6 


255.0 

6.30 

7 

12 

41.7 

12.2 

- 


68.4 

12.33 

- 


255.0 

12.4 

9 


These ingots were about 1 cm in diameter and 15 cm 
in length. The ingots were cut into small samples of 
approximately 2.5 cm in length. The ingots were cut into two 
halves along its axis and one half of it is used. Figure. 1 
shows a typical sample size. 


L 

1-icm 

T 



K- 2'S cm -H 



Figure 


1 . 
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Each sample is marked with an arrow which determines 
the face which, is to be polished for X - Ray exposure. 

The preparation of the sample for X-Ray exposure 
involves polishing and etching of the transverse 
surface (cross hatched in Fig. 1). Since the material is very 
soft, care has to be exercised at each stage of polishing. 
The initially cut sample is polishied with emery cloth. Then 
polishing is carried out in 1/0 , 2/0 and 3/0 emery papers. 
After fine emery polishing, the samples are polished in wheel 
polishers, fitted with a silk cloth, using Alumina as the 
polishing medium. This gives a fairly good, scratch free 
surface. 


The samples are etched and polished alternatively. 
At each stage, the microstructure is observed at a low 
magnification (50 - 100X). Two different etchants were used, 
one for 2% and 6% Bi aaloys and the other for 12% Bi alloy. 
The etchants used are as follows : 

% Bi Etchant 
2 and 6% KgCrgO^ 1 gm + Cone. HCl 6-7 ml 


+ Distilled water 86 ml. 
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12 % BeCi3.6H20 2 gm 

Cone. HCl 5 cc 

Distl. water 30 cc 
CgHgOH (95%) 60 cc 
(Taffs Reagent) 

The microstructure of the polished surface is 
observed at 50 - 100 X and the largest grain is located. 
Whithin this grain, the suitable spot for X-Ray exposure is 
located. For most of the sciinples, the grain at the geometric 
centre is selected, since alligning the X-Ray beam to fall on 
the geometric centre is much simpler. If a suitable grain is 
not available at the geometric centre, then an alternative 
grain is selected and its position is noted, by drawing the 
microstructure of the surface and measuring the approximate 
location of the spot. 

3.2 Obtaining the Laue back reflection patterns : 

The Laue back reflection film cassette has the 
dimension for a film of 8.3 X 10.1 cm. The X-Ray film is cut 
as closely to this size as possible and a centre hole is 
punched. The X-Ray film is enclosed in a black paper envelope 
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to prevent blackening due to light leaking through the sides 
of the cassette. A black paper is kept at the centre also, 
with a hole of the size of the collimeter to prevent any 
leakage of light. 

The samples were exposed to X-Rays at 20 mA and 30 
kV for a period of 7 to 8 hours. Cu target X-Ray tube is 
used. The intensity of the spots obtained is observed to be 
low and hence the current and voltage were increased to 30 mA 
and 40 kV for some of the samples. THe exposure time was also 
increased to 9 hours. 

3 . 3 Obtaining the input data : 


The input data are obtained from the LBR pattern by. 


(i) 

tracing the locations 

of 

the spots on 

a tracing paper 


to measure the X and 

Y 

coordinates 

of 

each spot , 

(ii) 

preparing a stereographic 

projection. 

to 

locate the 


zone axes and to measure the angles between them. 

The X-Ray film is kept over a glass plate and a 
tracing papar is fixed over the film. By Illuminating the 
plate from the bottom, the position of each spot is traced on 
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the paper. The X and Y axes of the film (assumed to be 
parallel to the edges of the film) are also drawn on the 
paper. The origin is plotted at the centre of the hole in the 
film. Thus a replica of the film is prepared. 

This replica is then superimposed over a graph paper 
with 1.0 mm division. The origin in the pattern is coincided 
with the origin on the graph paper. Then the X and Y 
coordinates of each spot is measuredwith an accuracy of +0.5 
mm using eye ball judgement. Appendix I shows the 
transparencies prepared for each film, with the X and Y 
coordinates of the spots. 

By using a Greninger’s chart, the r and 6 values 
of the spots are measured. These values are then transferred 
on a stereographic projection. THe spots belonging to one 
great circle are located and all possible great circles are 
traced, with each great circle containing atleast three 
spots. The great circles containing more than 4 or 5 spots 
are considered as low index zones. The spots lying at the 
intersection of more than 2 zones are considered as the 
prominent spots. The zone axis for each zone is located on 
the stereographic projection and the mutual angle between 
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these axes are measured. 

The proper selection of the desired number of spots is 
important for the indexing of the spots. The spots are 
numbered arbitrarily in the beginning. The spots lying at the 
intersection of two or three zones are renumbered first and 
the rest of the spots are numbered in that sequence. A 
minimum of six spots are selected such that they involve 
atleast two to three zones. The zone axes are labelled with 
alphabets . 
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CHAPTER 4 : COMPUTER PROGRAM 

4.1 Technique adopted : 

An iterative scheme of matching the iriterplanar 
angles between the planes causing the spots, obtained, 

(i) using their X and Y co-ordinates on the X-Ray film and 
(ii) with the assigned Miller indices, has been used in the 
Computer code developed. The scheme measures the angle 
between the spots using their X and Y co-ordinates, which is 
designated as , representing Experimental angles. Trial 

indices are assigned to the spots in a sequence and using 
these indices the mutual angles are measured using the 
crystallographic expression for the tetragonal crystal (eqn. 
8). These angles are designated as representing 

Theoretical angles. These two angles are compared such that 
they are within an allowable error limit (- 3.0 degrees). 

This scheme is carried out for all the spots. 

In addition to matching the interplanar angles between 
the spots, the conditions that, certain defined spots should 
belong to the same zone and the mutual angle of their zone 
axes should match with the given values are also carried out. 
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For each Laue film, a stereographic projection is prepared. 
From this the spots belonging to the Scime zone are identified 
by tracing the spots falling on one great circle. The zone 
axes are also located, on the stereographic projection. Then 
the mutual angle between these zone angles are measured on 
the stereographic projection, which in turn must match the 
ones calculated theoretically, using the appropriate 
expression in the program. 

The scheme also incorporates optimisation of the film to 
sample distance. This is done by minimising the sum of 
squared errors between the experimental angles and the 
theoretical cuigles (0^^- indices of each zone axis are 
also calculated and their mutual angles are calculated. 

This chapter explains the method of measuring the angles 
using the X and Y co-ordinates, selection of indices, 
indexing scheme, determining the crystallographic orientation, 
and adjustment of film to sample distance. 

4.1.1 Measurement of mutual angles : 

(i) Experimental angles : 

The experimental angles refer to the mutual angles 
between the spots obtained from their X and Y co-ordinates. 
The method adopted is the one developed by Christensen et al. 


17 
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[5,6]. The computation of angles associated with pairs of 
ref lectionsis accomplished in a perfectly straight forward 
manner. The geometry involved in the calculation of the 
angles is explained in Figure 2. The film plane is assumed to 
have the X and Y axes. The beam direction is taken as the Z 
axis. The specimen to film distance is given as D. Each spot 
falls on the film with corresponding X and Y values. The 

angle 24 . . is the angle between two reflections from the 

1 > J 

sample S. These two reflections have X and Y values as 
( X. Y, ) and ( X. Y.. ), with respect to the origin 0, which 
is the point through which the X-Ray beam passes through the 
film. From figure 2. a, 

9 9 9 9 

q = X.^ + Y^'^ + (1) 

2 2 2 ? 

= X^"^ + Y^"" + (2) 

L 3 = ( X. - Xj )^ + ( Y. - Yj (3) 

Figure 2.b gives the construction of the two spots’ locations 
and the sample, in which 



( 4 ) 
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L 


6 


2 




(b) 


~ ^ — L“ ^ ^ — ^L — ^ 

1 6 


~1 ^4 -1 ^4 

</’. , =: ]/2{ tan ^(— ,--) + tan ^ (— ^—-5— ) 

J i-i - jLi /-* 


( 6 ) 

(7) 


<^'’. . represents the angle between the spots, i and j. These 

aJ 

angles are measured for all the spots to be indexed. These 
are the experimental angles <^'^y i = 1.2, .... N and j - HI, 
i+2, . . . . , N, where N is the total number of spots to be 
indexed. 

(ii) Theoretical angles : 

For a set of indices assigned to two spots, the angle 
between them is measured using the crystallographic 
expression for the tetragonal crystal. These angles are 
designated as theoretical angles tp^. The expression for 
measuring the interplanar angle is 

(hih. + Vj)(c/a)2 + (Ijlj) 

COS ( 4'^^ ) = 

L(h^+ki)(c/a)^+1^3 [ (h^+k^) (C/a)^+1^3 


( 8 ) 
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where h,k,l are the indices assigned the spots i and j. 

(iii) Zone angles : 

The scheme incorporates angular match between the zone 
axes also. Hence for a set of indices (u.v.w. ) and (u.v.w.) 

1 3. iX X 0 

assigned to two zone axes, their mutual angle is given by the 
expression. 


(UiU. + v.vp/(C/a) + w.w. 


COS = 


[(u^+v^)/(C/a)^+ w^] t(u^+v^)/((^/a)2+ w^] 


( 9 ) 


4 . 1 . 2 . Selection of indices : 

The index search is restricted between ( 6 6 S ) to 
( 6 6 B ), since a judicious choice of spots on the film will 
normally ensure that their indices lie within this range. The 
total niimber of indices possibile in this range are 1331 
including ( 0 0 0 ) . For body centered crystal , from the 
structure factor expression, planes with (h+k+l)= odd 
do not produce any reflections. Hence all those indices with 
( h + k + 1 ) = odd are excluded from the range of indices. 

Considerable savings in search time can be achieved 
without any loss in generallity if the scalar multiple of 
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previously considered indices are excluded from the search. 
In general if(hkl)is considered, then ( nh nk nl ) need 
not be considered, where n is an integer. These conditions 
greatly reduce the number of index sets to be considered. 

The selected index sets within this specified limits are 
arranged from lower indices to higher indices, so that the 
lower indices are first assigned to the spots. For a given 
index set, ( h k 1 ), larger the d spacing, lower will be 
the indices. Based on this criteria, the index sets are 
arranged in descending order of d spacing, which is 
calculated for each index set from the expression. 






( 10 ) 


where h,k,l are the indices and a and C are the lattice 
constants . 


4.1.3. Scheme of indexing : 

The scheme finds out trial indices for the N spots, such 
that all the experimental and theoretical angles are within 
an allowable arror range. The error criterion is defined as. 


<P ^ - 4’ g 


< ERROR 


( 11 ) 
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The allowable error range is fixed as +3.0 degrees. All the 
spots are initial lised with the indices (011 ), which is 
the first index set to bo assigned. The heirarchy of indexing 
is explained in Figure 3. The scheme tests the angular match 
between a pair of spots at any stage. First it checks the 
angular match between spot numbers 2 and 1 . If the angle 
between these spots are not matching with the their 
experimental angles, then the indices of spot number 2 is 
changed. With this index, the angle between spot 2 and 1 is 
again checked. If a match is not found with all the possible 
index sets being assigned to spot number 2, then the indices 
of spot number 1 is changed to the next set in the list and 
the processing of assigning indices to spot 2 is repeated. 



Figure 3 

Angular match should be obtained between the spots 
in the tail end and head end of the arrows. 
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If a suilabie mcjtch is found between spots 1 and 2, then 
with their indices fixed, the scheme moves to the next spot, 
which is spot, number 3. Now the angular match between spot 
number 2 and 3 is checked, by assigning indices to spot 
number 3. With a pariticular index set assigned to spot 
number 3, if an angular match is obtained between spots 2 and 

3, then with that index set for spot number 3, the angular 
match will be checked between spot number 3 and 1. If the 
angular match between these two spots is not acceptable, then 
the indices of spot number 3 will be changed and the angular 
match between spots 2 and 3 will be checked. The indices 
assigned to spot number 3 will be changed, if the angular 
match between spots 3 and 1 are not acceptable, even though 
the match between spots 3 and 2 is satisfactory . If the entire 
range has been exhausted for spot number 3 and if still a 
match is not obtained, then the indices of the spot 
immediately before it will be changed, ie the indices of spot 
number 2, and the angular match proceeds from this spot 
onwards. Fixing the indices to a spot implies that the 
indices assigned to that spot satisfies the angular match 
with all the previous spots. For example assigning indices to 
spot number 5 would mean, angular match between spots 5 and 

4, 5 and 3, 5 and 2, and 5 and 1, This way the scheme 
proceeds the indexing of all the spots. 
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In addition to matching the interplanar angles, the 
scheme also assigns indices such that, they satisfy the zone 
condition, if a set of spots belong to one zone;. The numi)ers 
of the spots belonging one zone are obtained from the 
stereogaphic projection of the pattern. For eample in set of 
6 spots, if spots 1, 4, 5 and 2, 3, 6 are belonging to one 
zone, then the indices are assigned to the spots so that this 
condition is satisfied. In the indexing scheme, for eample if 
spots 1, 4, and 5 are belonging to one zone, then when 
assiging the indices to spot number 5, satisfying the angular 
match satisfactorily, with the indices of spot 1, 4, and 5, 
the condition, that they belong to the same zone is cheeked. 
This check is done by calculating the determinant of the 
matrix formed by the indices of the three spots. If the zone 
condition is not satisfied, then the indices of spot number 5 
will be varied and the angular match explained in the 
previous paragraph will be again carried out. This way the 
assigned indices satisfy the zone condition also for certain 
specified spots. 

The indices of each zone axis will also be calculated 
and the angle between these zone axes is matched with the 
given values. Essentially zone checking includes assigning 
the indices to the spots such that they are belonging to one 
zone and the angles between their zone zxes match with the 
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given values obtained from the stereographic projection. 

The sone checking condition is incorporated in the 
scheme to obtain a unique set of indices for the spots. The 
scheme explained, incorporates some of the features of manual 
indexing. The scheme fecilitates the possibility of obtaining 
a unique set, irrespective of the sequence in which the spots 
are numbered and fed. Since it is quite possible to obtain a 
completely different set, which also satisfies all the 
required conditions, when the sequence in which the spots are 
numbered is changed, the zone checking condition is 
important. It is observed, that inclusion of this condition 
results in a unique set of indices for the spots, 

irrespective of the sequence of numbering them. 

The scheme explained is one of trial and error method 
and the proper selection ofspots for feeding into the program 
is very important for obtaining a set of indices, satisfying 
all the required conditions. The spots which lie at the 
intersection of two or three zones will be lablled first. The 
selected spots will include a miniminn of 2 zones, so that 
from the indices of these zone zxes, the orientation of the 
crystal can be determined. 

4.1.4. Determining the crystallographic orientation : 

The scheme, after indexing the required number of spots 
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gives 

the 

indices 

of 

the 

zone 

axes 

also. 

From 

these 

indices 

of 

zone axes , 

the 

orientation 

can be 

determined 

easily . 

The 

samples 

used in 

the 

present 

work 

were 

Sn-Bi 


alloys, which have their primary dendrites roughly parallel 
to the growth direction, which is [ 1 1 0 ], as reported in 
the literature [1,2,3]. The angle between the indexed zone 
axes and [ 1 10] is calculated . The method of determining 
the orientation is as follows: the stereographl c projection 
of the pattern is fixed on a Wulff net, such that the X and Y 
axes of the projection coincide with the E-W and N-S 
directions of the net respectively. The pole of the axes is 
first rotated with respect to the N-S axis to bring the pole 
on the curcumferance of the net. (refer to Figure 4). The 
angle by which the pole is rotated about the N-S axis is 
noted. The pole is then made to coincide on either the North 
or South pole of the net. With this position of the pole, the 
small circle which is at the required angular distance from 
the pole is traced. Few points are marked on this small 
circle. Now the pole is rotated back with respect to the axis 
perpendicukar to the projection. The pole is again brought to 
its initial position by rotating with respect to the E-W axis 
by the same nxamber of degrees, by which it was rotated 
before. The points which are marked on the small circle are 
also rotated by the same number of degrees about the N-S 
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axis. The arc containing the new position of these points is 
traced, which inturn will be at the requird number of 
degrees from the zone axis. This process is repeated for the 
next zone axis also. The intersection of the two small 
circles traced for each pole, fixes the location of the 
desired zone axis ([110]). Now the angle between this 
axis and the centre is measured by alligning the two poles to 
lie along the equator. In terms of this angle, the 
orientation of the crystal is expressed. 

4.1.5. Optimising the distance : 

The best estimate of the film to sample distaince, D is 
determined by minimising the sum of suqares of the errors 
between the experimental angles 4 ^^ and the theoretical 
angles 4 ^^, obtained with the final indices of the spots. The 
method adopted is the one described by Christensen et al. 
[5,6] and is as follows : 

The mean square error, Z, is defined as, 

Z(D) E ( 4 > p.(D) - 4 > ^ (12) 
i=l j-i+i ^ ^ 

where n is the number of spots indexed. The error is a 
function of the distance D. This error function is 
approximated by the quadratic. 
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Z(D) - 80 -> a^D + 


(13) 


The error is minimised wilh respecl lo D when, 


" a. 


D = D = 


2 3 , 


if 28^ > 0. 


(14) 


The co-efficients a , a, , and a„ are determined by a 3 point 

0 1 ii^ 

fit of the equations of the form : 


Z. = a~ + a*D. + a,,D.‘ 
1 0 1 1 2 1 


(15) 


where i = 1,2, or 3 and is obtained by letting D = D^. The 
initial values of D^, Dg and are taken as D, D - 0.1 and 
D + 0.1 cm respectively. With the stcindard value of 3 cm, 
this scheme provides an initial variation in D of +1 mm. 
Substitution of the co-efficients 3 ^, a^ and 3 ^ in equation 
15, gives. 


D* = 


ZiCD^-Dg) + ZgCDg-Dj) 4 Z3(D,-D2) 


(16) 


With the initial value of D as 3.0 cm, a new value of D will 
be calculated using the expression 16. If these two values 

3jC 

are within 0.1 mm of each other, then D is taken as the 
optimised distance. Else, with D as D, expression 16 is used 
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1,0 obtain a new value of D . This iterative procedure is 
continued until the two successive calculations are withim 
0.1 nun of each other. The best value of D obtained is used to 
calculate the optimized experimental angles. 

The optimised distance is expected to be within 3.2 to 
2.8 cm. If a value of best distance within this range is not 
obtained then, the scheme calculates the distance, within 3.2 
to 2.8 cm which has got the least mean square error . The 
reason for putting such a restriction is explained in the 
next section. 

4.1.6. Criteria for fixing the error limit : 

The computer code is developed with an allowable error 
limit of + 3.0 degrees, in matching the interplanar angles. 
The error limit is set based on the errors likely to have 
been introduced during the measurement of input data. These 
are as follows : 

(i) The spots on the X-Ray film generally do not have an 
exact location. The identification of the eaxct location is 
difficult, if the spots have fuzziness. 

(ii) The experimental angles are calculated with the 
standard value of D as 3.0 cm. But, while fixing the sample 
for X-Ray exposure, the distance of seperation between the 
sample and the X-Ray film could be any where between 2.8 to 
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3.2 cm. With sufficient care it will be possible to fix the 
distance within this range. Hence this limit is kept in 
optimising the distance of seperatiori. But since, the 
experimental angles are calculated initially with D as 3.0 
cm, there will a certain amount of error introduced in the 
experimental angles. 

(iii) The X-Ray beam may not be passing through the centre 
of the hole punched in the film and hence the X and Y axes 
for the spots on the film may not be parallel to the edges of 
the film. 

These errors can vary the measured experimental angles. 
The effect of X and Y values on the experimental angles is 
calculated. A set of 5 X and Y values were selected and they 
are varied 1.5 mm on both sides of the selected value. With a 
fixed value of D, the experimental angles are calculated for 
each value of X and Y. These angles are given in the Computer 
output No.l, in Appendix III. A maximum error of about 1.0 
degree is observed between the extreme ends of the X and Y 
values. 

The effect of variation in experimental angles due to 
variation in D value is also calculated. For the same values 
of X and Y, experimental angles are calculated with the 
distance being increased in steps on 1 mm. These angles are 
given in computer output No. 2 in Appendix III. Here a maximum 
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error of about 3.2b degrees is observed. 

The error due to the X-Ray beam not passing through the 
centre of the punched hole is eliminated by fixing lead 
markers in the X-Ray film cassette. This fixes the X and Y 
axes of the film without ambiguity. This method has been 
tried for sample number 6.4 M. 

Based on these calculation, the upper limit for the 
allowable error is kept at 3.0 degrees. 

4.2. Structure of the Program : 

The computer code developed has been divided into number 
a of subroutines, each performing a specific task. The 
subroutines are arranged in the sequence, in which the tasks 
are to be performed. They are called from the main program in 
the same sequence. This chapter explains the function of each 
subroutine. The input data to be fed and their sequence is 
also explained. Flow cahrts are givan for each suroutine in 
Appendix IV . 

4.2.1. Input Data : 

The input data required are listed below in the same 
order in which they are fed. 

1. Total number of spots to be indexed (N) 
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2. X and Y co-ordinates of each spot obtained from the X-Ray 
film, in cm 

3. Total number of zones to be checked (NZ) 

4. Numbers of spots belonging to each zone 

5. Nutual angle between each zone. 

Additional input data required are, error allowed in the 
angular match (ERROR), the maximum and minimum limits of the 
Miller indices to be assigned to the spots for index search 
(MAX and MIN), C/a ratio of the crystal (RATIO), the distance 
of separation between the film and the sample (BIST). These 
data are given in the program itself, since they remain 
unchanged for all the samples. 

A maximum of 15 spots can be indexed. This limitation is 
not absolute, but can be changed as required, by redefining 
the storage space, for the appropriate variables. The X and Y 
co-ordinates of each spot to be indexed should be obtained 
with an accuracy of +1.0 mm. The procedure for obtaining the 
X and Y values has been explained in section 3.3. A stereogram 
of the pattern is to be prepared to identify the spots 
belonging to one zone. All the possible zones are identified 
and the zones which contain more than three spots are 
considered as important zones. Each zone is located on the 
stereogram and their mutual angles are measured. The spots are 
selected and numbered accordingly as explained in section 3.3. 
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The set of indices obtained is expected to be Scune for any 
order in which the spots are fed. But the number of iterations 
in each case may be different. For sample number 12.4, the 
spots are fed in different sequences and it is observed that 
the set of indices obtained remain the same, irrespective of 
the sequence. 

4.2.1. Subroutine TRIANG( lA, IB) : 

This subroutine calculates the angle between any two 
planes using their indices with the help of the 
crystallographic expression, for tetragonal crystals (eqn. 8). 
The counters lA and IB represent the spot numbers whose mutual 
angle is required. The latest indices assigned to the spots 
are used for calculating the angles. The C/a ratio value is 
transferred from the maitx program through COMMON statement. 
This routine is called from subroutine TEST, while the angxilar 
match is being carried out. The calculated angle is returned 
as TANGLE{IA, IB) representing the trial angle between the spot 
number lA and IB. When all the spots are indexed, this routine 
is called again from the main program and the mutual angle 
between the spots are measured and stored as theoretical 
angles. THis routine is not called from any other subroutine. 

4.2.3. Subroutine EXPANG ( D ) : 

This subroutine calculates the mutual angle between the 
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spots using thoir X and Y co-ordinates. The steps involved in 
calculating the angles are explained in section 4.1.1. The 
quantity D represents, the distance of seperation between tlie 
sample and the film. The angles can be measured for any given 
value of distance. This subroutine is called in the beginning 
from the main program with D as 3.0 cm and the calculated 
angles are returned as EANGLE{I,J). These angles are stored as 
EANG(I,d) permanently, representing the experimental angles. 
THis subroutine is also called from subroutine ERRO, since the 
mutual angles are required for different values of distance D, 
to calculate the mean square erron bwteen the experimental and 
theoretical angles, during optimising the distance of 
seperation. 

4.2.4. Subroutine ZANGL(MZN, IZ , ZANG) : 

This subroutine calculates the angle between the zone 
axes, using the appropriate expression for the tetragonal 
crystal, (eqn. 9). This routine is called from subroutine 
ZONEM for matching the angles between the zone axes. The 
calculated angle is returned as ZANG, rei resenting zone angle. 
The counters MZN, and IZ represent the nimbers of the zones 
whose mutual angle is required. The indices of the zone axes 
are transferred through COMMON statements from the calling 
subroutine. The C/a ratio is transferred from the main 
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proeraiD. After the indices of al] the spots are fixed, the 
zone angles are calculate using the final indices of the zone 
axes. This is carried out from the? subroutine ZONE. 

4.2.5. Subroutine SELECT : 

The indices to be assigned for the spots is restricted 
between (555) to (555). Hence all the possible indices 
which could be assigned to the spots are first selected within 
this range and they are stored permanently as 1H,IK,IL. The 
selection of such a list is carried out in this subroutine. 
The selected indices are then rearranged from lower indices to 
higher indices based on the value of their d spacing. The 
starting index set is ( 000 ), which obviously will not be 
assigned to any spot. Three counters H,K,L are used in 
selecting the indices. The counter is first variedfrom 0 to 
the maximum limit in steps of 1. Whem L reaches the maximunt 
limit, its value is changed to -1 and subsequently it is 
decreased in steps 1 upto to the minimum limit. After the 
counter L has been changed from +5 to -5, . its value is 
reassigned as zero. This represents that the counter L has 
undergone one full cycle. Now the counter K is increased from 
zero to +1. With this value of K, the counter L again under 
goes one full cycle. Then K is increased in steps of 1 upto 
the maximum limit and for each change in K, L undergoes one 
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full cycle. When K completes one full cycle, ( ~b to +5), the 
counter H is varied from -5 to +5, which is the last counter 
to b<; varied. With every change in counter H, both K and I. 
have t.o undergo their respective cycles, to exhaust all 
possibile combinations of the Miller indices. Counter H needs 
to undergo only one cycle. When the values of H,K, and L 
reaches the minimum limit, it represents the end of the 
selection and all the counters are reassigned as zero. The 
total number of indices possible within this range are 
11x11x11- 3 = 1330. Certain index sets are eliminated in 
this range, as they either represent the same planes or are 
otherwise not permissible. Eliminations of such indices is 
carried out in the routine LAOE and the eliminations are 
explained in section 4.2.6. Each index setbefore being 
selected to be in the list is tested for these elimination 
conditions. After all the eliminations, the list contains a 
total of 578 index sets for (555)to(555). The index 
set ( 000) is also included in the list, as the last entry. 

The selected indices are rearranged from lower indices to 
higher indices. The larger the d spacing for a given index 
set, the lower will be the indices. On this basis, the d 
spacing for all the selected index sets are calculated and the 
index sets are rearranged in descending order of d spacing. 

The indices from this list are assigned to the spots one 
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by one. The sfdected index sets are made available to the main 
program and the required subroutines. 

4.2.6. Subroutine LAUE : 

Each index set (to mean one combination of H,K and L), 
selected by the routine SELECT, is transferred to this routine 
to check, whether it is to be rejected by any of the 
elimination conditions. The elimination conditions checked are 
(i) the sum ( H + K + L ), should be an even number. 

The structure factor expression is, 

F = E f (17) 

( = 1 

where f is the atomic scattering factor, h,k, l are the 
Miller indices of the plane anf u,v. n' are the atom positions 
in the unit cell of the crystal. For a body centered crystal 
the atom positions are 0,0,0 and 1/2, 1/2, 1/2. Hence the 
structure factor expression for a body centered crystal is. 


F = 

f( 1 + e-^'" 

h +k +i :> ^ 

(18) 

for 

( h + k + 1 

) = even, F = 2f ,I o F^ = 4f^ 

(19) 

for 

{ h + k + 1 

) = odd, F = 0 and 1 = 0 

(20) 


This implis that, there will be no reflection from planes 
with h h I - odd. This is true for all boby centered 
crystals irrespective of the system to which they belong. 
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(li) reduction of the Indices to be tried can be achieved by 
noting that same value of angles result from ali scalar 
muJtipies of an index set. Thus, if ( 1 1 0 ) is considered, 
then (220) or (330) need not be considered. In general 
if (hkl ) is considered, then ( nh nk nl ) need not be 
considered, where n is an integer. 

These conditions are selected for each index set selected 
and the decision is transferred through the logical variable 
EXTIN, to subroutine SELECT. If EXTIN is true, the particular 
set will not be included in the list and the selection 
procedure will continue to the next index set. This routine is 
called only from subroutine SELECT and not from any other 
routine. 

4.2.7. Subroutine INDEX : 

This subroutine does the ’Indexing’ of the spots. Each 
spot is assigned with an index set which satisfies the 
interplanar angular relationship as well as zone relations. 
The spots are assigned indices from the list selected by the 
routine SELECT. The routine TEST is used for checking the 
angular match. Routine INDEX is called from the main program 
only once. All the input values are transferred to it through 
comoion statements. 

The important variables used in this routine are. 



41 


(i) counter ]; this represents the number of the spot to 
which the indices are assigned at any stage, 

(ii) Logical variable ACCEPT; this takes the value true, if 
the indices assigned to the spot satisfy the angular match 
with all the other previously indexed spots. This variable is 
obtained from the routine TEST, since the angular match is 
checked in that routine, 

(iii) Logical variable RANGE; if all the indices in the list 
have been tried for a particular spot, then the variable 
RANGE takes the value false. Thus, all the indices have been 
tried out for that spot and since no match is obtained, the 
indices of the previous spot has to be changed and the 
sequence of trying indices for the spot is repeated . 

All spots are initially assigned the indices (011 ) . 
The counter I decides at any stage which spot is being 
indexed. This subroutine has the provision, to see the 
intermediate values of the indices assigned upto a certain 
instant. Execution is paused after each 50000 iterations 
(this number can be changed), and the routine displays on the 
screen, the indices assigned to the spots upto this stage and 
asks if further execution is desired. If number 1 is keyed 
in, then the execution is continued and for any other numeral 
the execution is stopped and it comes to an abrupt end. 
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4.2.8. Subroutine TEST(J) : 

This subroutine checks whether the angle between any two 
spots is within the allowable error limit ( j^y.0 degrees ). 
The counter J^l represents the number of the spot which is 
undergoing the checking. The check is carried out using the 
routine TRIANG ( see section 4.2.2 ), which calculates the 
interplanar angle with the indices assigned at any stage. 
This angle is then compared with the experimental angle 
between the corresponding spots, measured in the beginning 
from routine EXPANG. If the difference in the angles is 
within the allowable error limit, the indices assigned to the 
spots are accepted. The logical variable ACCEPT is used for 
this. Routine TEST also checks whether a given spot requires 
zone cheking and if so it decides whether the indices satisfy 
the zone condition. The subroutine ZONECK is called for 
checking the zone condition. 

The indices assigned to a spot should satisfy the 
angular match with all the previous spots. The counters MN 
and MM are used for this check. The counter MN initially 
takes the value of J and MM is J+1. The value of MM is kept 
at J+1. Now the match between spotn\aibers MN and MM is 
checked. If it is acceptable, MN is decreased by 1 (thus MN 
becomes J-1) and the angular- match is checked between the 
current value of MN and MM. In this way, counter MN is 
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decreased upto 1. The indices assigned to spot number MM 
should satisfy the angular match with all the previous spots 
represented by the counter MN. Then only the logical variable 
ACCEPT takes the value true. If the match is not satisfactory 
at any stage then ACCEPT becomes false and the indices 
assigned to spot number MM should be changed. 

Zone checking will be carried out only if ACCEPT is 
true. Zone checking is done only for specified spots already 
identified in the beginning. These spots are designated as 
zone points. When accept becomes true, the spot nuBiber MM is 
checked whether it is one of the zone points. This condition 
is determined by the logical variable ZONE. If ZONE is true, 
then the subroutine ZONECK will be called, in which the zone 
checking is carried out. If the zone condition is satisfied, 
the variable ZONIC becomes true. Depending upon the value of 
ZONEC, the variable ACCEPT finally takes the value true or 
false. If ZONEC is true, ACCEPT is also made true. If the 
final value of ACCEPT is true, it represents that the indices 
assigned to the spot MM (J+l), satisfies the angular match 
with all the previous spots and also the zone condition, if 
it is a zone point. 

4.2.9. Subroutine ZONECK (M) : 

This subroutine has to check two conditions, 
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(j) does the spot belong to the same zone as specified ? 

(ii) if so, do the calculated values of the angle between 
various zone axes, match with the given values, within an 
error of + 3.0 degrees. 

Condition (i) is checked in this routine itself and 
subroutine ZONEM is called for condition (ii). First the 
number zones to which the spot is common is determined. This 
is decided by the variable NZSB. Then for each zone, the 
condition that they should belong to the same zone is 
checked. If a given spot belongs to more than one zone, the 
condition (i) should be satisfied for all these zones. The 
decision is determined by the logical variable ZONEC. 

Once condition (i) is satisfied, the calculated angles 
between zone axes should match with the given values obtained 
from the stereographic projection. This checking is required 
only if the spot belongs to zone numbers 2 or higher zone 
numbers. Hence, first the maximum number of the zone to which 
the spot belongs is determined. This is given by the variable 
MZN. Since angular match for zone axes is carried out in 
subroutine ZONEM, for MZN equal to 1, this subroutine will 
not be called. Otherwise, the zone axis of zone MZN is 
calculated first. The indices of zone axex are designated as 
Ml ( I ) , M2 ( I ) , M3 ( 1 ) , where I is the number of the zone . The 
counter IZ is used to number the other zone which is checked 
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with zone MZN. IZ initially takes the value MZN-1, and 
routine ZONEM is called. The decision is made by the logical 
variable ZONMAT. If ZONMAT is true, then IZ is reduced by 1, 
and the angukar match between zone IZ MZN is carried out. In 
this way IZ is reduced upto zone number 1. If at any stage 
ZONMAT becomes false, further checking is not carried out. 

Depending upon the value of ZONMAT, the logical variable 
ZONEC is reassigned the value true or false. If ZONMAT is 
true, ZONEC is assigned the value true. ZONEC is then 
transferred to the calling routine TEST. 

4.2.10. Subroutine ZONEM{IZ) 

\ 

The counter IZ has the same meaning as in the routine 
ZONECK. This routine ZONEM, first calculates the zone axis of 
the zone IZ. Then with these indices, the angular match 


between 

the 

zone IZ 

and MZN 

is checked . 

This is 

done 

by 

calling 

the 

routine 

ZANGL, 

which gives 

the value 

of 

the 


calculated zone angle as ZANG. This value is then compared 
with the given value for those particular zones. If the match 
is acceptable, ZONMAT becomes true. This routine is not 
called from any other subroutine. 


4.2.11. Subroutine ADJOST{DI) : 

This subroutine optimises the distance of separation 



46 


be'Lween "the sample and the X-Ray film. The initja] distance 
is assumed to be 3.0 cm. This distance is optimised by 
minimising the difference between the theroreticai and 
experimental angles. 

The scheme of optimising the distance is explained in 
section 4.1.4. The counter DI represents the initial value of 
the distance. The mean square error for this distance is 
initially calculated. The calculated error is stored as ZE. 
The distance is varied in steps of Imm on both sides and the 
new value of distance is calculated. If the calculated value 
is within 0.1 mm of the previous value, this new value is 
taken as the optimised distance. Else, the original value is 
changed to the new value of distance and the process is 
repeated . 

A limit is set, within which the optimised distance is 
expected to be found. This limit is 2.8 cm to 3.2 cm. The 
reason for fixing this limit is explained in section 4.1.6. 
If the above scheme is not able to find the optimised 
distance within this range, then the program will select a 
distance within this range which has the least mean square 
error . 

The logical variable FINAL is used to detennine whether 
the distance range is exhausted or not. The variables DSTAR 
and ZERROR thakes the value of the optimised distance and the 
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error associated with that distance, respectively. 

4.2.12. Subroutine ERR0(DK,1)Z) : 

This subroutine calculates the mean square error between 
the theoretical and experimental angles for a given value of 
distance using equation 12. The variables DK and ZK represent 
the value of distance and the mean square error respectively. 

The value of experimental and theoretical angles are 
transferred from the calling routine ADJUST. 

4.2.13. Subroutine ZONE : 

After all the spots are indexed, the zone axes of each 
zone and the mutual angles of zone axes are calculated in 
this routine. This routine is called from the main program. 

The zone axis is calculated using the indices of any two 
spots belonging to that zone. The indices of the zone axes 
are M1(I), M2(I), M3(I), where 1 is the number of the zone. 
The indices Ml, M2, and M3 are given as, 

Ml = C k l . - I * h (21) 

^ I J i 1 

M2 I h . - h .X I ( 22 ) 

I j i j 

M3 = C h k . - M * h (23) 

1 i J 

where and and ij are the indices of spots 
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I and j . The zone axen are calcula'ted for all the zones. 

The mutual angle between the zone axes are calculated by 
calling the subroutine ZANGb. The calculated angles are 
stored as ZANGLE(I,J) where I and J represent the zone 
numbers . 

4.2.14. Main Program : 

The main program first reads all "the input values fed in 
the specified order, from the file IN. IN. From the numbers of 
the spot belonging to one zone, the zone points are 
identified for each zone. Zone point is the number of the 
spot belonging to that zone, which is selected for carrying 
out the zone check. ITor each zone, three spot numbers are 
given. The maximum of these three numbers is taken as the 
zone point, for this zone. For example, if spots 2, 3 and 5 
are belonging to one zone, spot 5 is taken as the zone point 
for that zone, because zone checking need not be carried out 
for spots 2 and 3. The counter ZP(1), represents the zone 
point for zone I. 

After identifying the zone points, the zones are 
renumbered as 1,2,3 etc according to the zone points which 
are arranged in ascending order. The logical variables ACCEPT 
and RANGE, to be used later in the other subroutines, are 
given the initial valuta false and true respectively. All the 



49 


spots are gjven the intial indices (011 ). 

The subroutine EXFANG(3,J) is called with the intial 
value of distance (3.0 cm) and the experimental angles are 
stored as EANG(1,J) with 1 and J representing the spot 
numbers . 

The main program calls the various routines in the order 
in which their tasks are to be performed. The routine INDEX 
is called first, since the indices of the spots are to be 
fixed first. INDEX transferrs the indices of the all the 
spots satisfying the required conditions. With these values 
the theoretical angles are calculated and the other 
subroutines ADJUST and ZONE are called. The optimised 
distance DSTAE is used finallly to calculate the 
experimental angles. 

The subroutine PRINT is called finally which prints both 
the input and output values in a structured format. 



50 


CJIAPTKU 5. WiSULTS AND DISCUSSIONS 


Usirifi the computer code developed, eight Laue back 
refJert.jori patteriir. for Sn-Bi alloys were Indexed and the 
orientation of each sample was determined, from the Indices 
of the zone axes. The location of [1 10] direction is 

located on the stereographic projection for each sample and 
the angle between [110] and the centre is measured. The 
set of indices obtained for each sample from the computer 
program are given in Appendix III. 

Table. 2 gives the results obtained from the computer 
program for each sample. 

TABLE 2. RESULTS 


Sample 

No. 

Computer 

output 

No. 

Number of 
spots 
indexed 

Number of 
zones 
indexed 

Maximum 
error in 
angles (deg) 

2.1 

3 

5 

2 

2.5859 

3.1 

4 

5 

2 

1.8539 

4.3 

5 

5 

2 

2.2688 

6.2 

6 

6 

3 

2.2354 

6.30 

7 

6 

2 

3.9867 

6.4 

10 

6 

2 

3.1596 

6.4 M 

8 

5 

2 

0.8112 

12.4 

9 

7 

3 

1.8656 
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Table 3. gives the orientation of each sample calculated 
from the stereographic projection of the pattern. 


TABLE 3. ORIENTATION OF THE SAMPLES. 


Sample 

No. 

Indices 
of zone 

Angle between [110] 
and zone axes (deg) 

1 

Angle with' 
[110] and 
X-Ray 
beam (deg) 


CZ A 

measured 

calculated 

2.1 

[ill] 

89.0 

90.00 

54.0 


[311] 

63.0 

63.83 


3.1 

[10 3] 

112.0 

111 .59 

82.0 


[110] 

0.0 

0.00 


4.3 

[ 0 1 1 ] 

51 .0 

51.61 

84.0 


[13 3] 

38.0 

37.39 


6.2 

[010] 

134.0 

134.95 

49.0 


[371] 

111.0 

111.70 



fill] 

90.0 

90.00 


6.30 

[111] 

89.0 

90.00 

76.0 


[111] 

158.0 

158,84 


6.4 M 

[111] 

22.0 

21.08 

27.0 


[011] 

51.0 

51 . 67 


6.4 

[16 3] 

122.5 

121.83 

not 


till] 

159.5 

158.84 

known 

12.4 

[12 0] 

70.5 

71.54 

88.0 


[110] 

0.0 

0.00 
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zone axer. , indexed by the program and | 1 10]. The method 
of locating the f 1 3 0 J axis on the projection is explained 
i n secti on 4.3.3. 

The published [3,2,3] growth direction for pure tin and 
tin-rich bismuth alloys is [ 330]. Hence, this axis is 
located on the stereographic projection and the angle between 
this axsi and the centre of the projection is measured. It is 
observed from the results, that the [310 ] direction is 
not coinciding with the centre and the angular difference 
between [330] and centre is observed to be between 27.0 
degrees to 88.0 degrees. The least angular difference is 
obtained for the sample 6.4 M (computer output no. 8, 
Appendix III ). The X-Ray pattern of this sample has been 
obtained, after fixing lead markers in the inner cover of the 
film cassette. The maximum angular mismatch, for this sample 
is also found to be the least (0.81 degress). 

In the light of these results obtained after fixing the 
lead markers in the film cassette, it may be worth while to 
take the LBR patterns for all the other experimental samples 
again. 

It is clear that, the growth direction for the primary 
dendrites for the experimental samples of tin-rich alloys has 
not been found to be [ 3 1 0 3 as published in literature. 
Tlie reason for this difference is not obvious at the moment. 
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Before conciudiriff anything about the growth direction for the 
experimental samples, a series of confinnatory tests need to 
be taken. Sucij a.s, 

(i) the scheme for indexing of the LBR patterns using 
computer should be modified, so that the zone axes are 
directly indexed rather than the spots. 

(ii) another set of LBR patterns shuold be taken, in using 
the film cassette with the lead markers 

(iii) specific grains need to be selected for LBR patterns 
atleast in some samples, which have sections in both the 
transverse and the longitudinal cut surfaces. This will 
permit the experimental measurement of the primary dendrite 
and its angle with thegrowth direction of the ingot, using 
optical microscope and micrographs. 

Finally, as an independant confirmation for the 
crystallographic orientation of the samples, electron 
diffraction patterns can be taken using Transmission 
Electron Microscope. 



CHAPTER 6. CONCLUSIONS 


(i) A computer program in FORTRAN 77 has been developed 
for indexing Laue back reflection patterns for body centered 
tetragonal crystals. The program indexes the spots of the 
Laue pattern and tiieir zone axes. 

(ii) The computer code developed has been used for indexing 
eight Laue patterns obtained from directionally solidified 
Sn-Bi alloys. The crystallographic orientation of each sample 
has been determined. 

(iii) The primary grovrth dirction for Sn and Sn-rich alloys 
which has been reported as [ 1 1 0 ] » has not been confirmed 
for the experimental samples used in the study . 
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CHAPTEfi Y. GCOPB; P“0K FURTHER WORK 

The present work carried out can be modified or extended in 
two stages, 

{ i ) in the experimental work of preparing the samp] es and 
obtaining the Laue back rection patterns, 

(ii) modification of the computer scheme employed. 

7.1. Modifications in the experimienta] work. 

The samples require polishing and etching for X-Ray 

exposure. The tin-rich bismuth alloys used were at present 
polished in wjhee] polishers using alumina as the polishing 

medium. Since these alloys are very soft, mechanica] 

polisihng in some cases results in a distorted layer 
affecting the structure of the surface, which ultimately 
results in a ring pattern after X-Ray exposure. 

Alternatively, continuous polishing and etching using 
che-rnoi'S leather has been tried and marked improvement in the 
quality of polishing is not observed. Electropolishing of the 
sample can be employed, which removesthe top layer without 


5b 
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affectinfi the structure. Since, a suitable set up was not 
available in the laboratory and within the available time 
frani'' it, was not possiijJe to fabricate one, this method is 
not tried. 

The Laue film cassette requires markers for locating the 
X and Y axes on the film without ambiguity. Lead markers are 
provided in the inner cover of the cassette and an X-Ray 
pattern is obtained for sample no. 6.4 M only. This fixes the 
origin on the film clearly and hence the X and Y co-ordinates 
can be measured with considerable accuracy. This improved set 
up may be tried for other samples as well. 

The X and Y co-ordinates of the spots were read using an 
ordinary graph paper with an accuracy of 1mm. A set up can be 
fabricated which has provisions for illuminating the pattern 
from the bottom andfor direct measurement of the X and Y 
co-ordinates. A set up similar to the one used for reading 
the lines of Debye-Scherrer patterns can be employed. 

These modifications will considerably increase the 
confidence with which the input data are obtained, which 
ultimately will reduce error in angular mismatch and also the 

computer time. 
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Y.2 Mod j f i oali onr. of the computer schemer : 

The present scheme first fixes he indices of the spots. 
Ijijiurd on the indices of the spots, th<r indices of the zori<r 
axes are fixed. For determining the crystallographic 
orientation of the sample it will be enough to find out the 
indices of the zone axes alone. Hence, alternatively direct 
indexing of the zone axes could be employed. 

The alternate scheme which could be employed is as 
follows : 

From the stereographic projection, one or two spots 
which are lying at the intersection of atleast three zones 
shuold be indent if ied. Then three to four zone axes should 
also be identified with their mutual angles being measured. 
The input data will be the X and Y co-ordinates of the spots 
(if two spots are used), the angle between the zone axes. The 
scheme will assign indices to the spot(s) (which is expected 
to be of low indices). Then simultaneously the indices of the 
zones will also be fixed. The angular check will now be 
between a plane and a direction. The indices of the normal to 
a given plane for the tetragonal crystal can be obtained as 


follows, 


h 

u 


k 

V 




( 23 ) 
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With thfi Indicefi of the plane being converted as the indices 
of a direction, the angle between the zone axes and the 
planes canbc measured, which will checked with the given 
values from the stereographic projection. 

This scheme will require two lists of indices, one for 
the planes (spots) and the other for the directions (zones). 
The selection procedure used in the present scheme Itself 
can be employed, for preparing these lists. 

Modification of the current scheme to the above 
mentioned one is not carried out in this project due the time 
limitations. It will be worth while to try this scheme, since 
the indexing of the patterns by this scheme is expected to 
give interesting results. 
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APPEFDIX. Ill 


COMIUTM OUTPUTS 
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14. 1953E 

IE. 65600 

1 .93594 

1 .53931 


.argest difference between the Angles ■ 3,9867 
imallesT difference between the Angles * ,£813 


ipots 

1 

£ and 

3 belong to tone 



c 7 1 7 3 


(pots 

4 

5 and 

1 belong to tone 

c 7 7 1 3 






Theoritical 

Exprimental 

^ngle 

between C 

Till and C 7 T 

1 3 » 97.407150 

95.500000 


optimised distance ■ £,8£16 

lean square error in angles with Optimised Distance ■ 45,88&£84 


OUTPUT NO. 7 




aapl« Number: fi.4H 


10. 7.1990 


Number of points indexed » 5 


Coordinates of the Points 


(X> 

(Yl 

.10 

.95 

1 .95 

0 

cu 

3.05 

.55 

■2.25 

1 .20 

■1.35 

-.05 


The index search is between ( £ 6 6) and ( £ £ £} 


Indices for Spot^ 1 
Indices for Spot 2 
Indices for Spot 3 
Indices for Spot 4 
Indices for Spot 5 


<011} 
C T 3 a > 

{ 4 5 5 1 
(233) 
(£35) 



Error range allowed “ 3.00 

deg. 



Points 

Expt . Ang Thert . Ang Opts . ang 

a 3 4 

Diffrnce. 

(3-2) 

Opt .Oiff 
(4-3) 

» 



16.10S7S 
ai .8ia34 
1 6 . 4354E 
15.88153 
14.450SS 
3t .68161 
31.91188 
. . 48 . B43B3 
■ 34 . 9833 ^.' 


15.94569 
aO . 9537a 
17.69884 
15.66975 
13.64668 
31 .43349 
31.61542 
38. 65854 
34'-2S3^S 

if ^ Li iilT^ ' r V , 


15. 
81 . 
17. 
15. 
14. 
31 . 
31 . 
39^ 
■■34''. 


81857 
33066 
98196 
31705 
187J59 
19319 
1 1368 
30881 


,16109 

.85868 

.73657 

.15178 

.80364 

.44818 

,81.64# 



it. 

:pm 

■ '5%%* 


.13318 
.37695 
.8631 1 
w'3587f 

.54088 

.'fe405|. 

.58:47#,: 






3««llest diffftrenc« betw«en the Angles » .1331 


Spots 

1 

2 and 

5 belong to zone 

C 1 1 T 3 


Spots 

1 

4 and 

3 belong to zone 

C 0 1 T 3 






Theorit ical 

Exprinental 

Vngle 

between t 

1 1 7 3 and C 0 1 

T 3 - 41 .262600 

39.060000 


Optimised distance • 3.1112 

lean square error in angles with Optimised Distance » 2.083834 


OUTPUT NO. 8 



8«Rpl« Number: IS. 4 


5. 7.1990 


Number of points indexed » 7 


Coordinates of the Points 
CX) (Y) 


.35 

.90 

1.75 

.05 

“.05 

“1.15 

-1 .05 

1 .95 

.60 

£.35 

t .90 

1 .30 

“1 .£5 

.68 


The index search is between { 5 5 5} and (555) 


Indices for Spot 1 ( 0 0 £ } 
Indices for Spot .2 * t £ 1 S ) 
Indices for Spot 3 ■ { 1 1 2 ) 
Indices for Spot 4 « ( 2 1 S ) 
Indices for Spot 5 ■ ( t 1 4 ) 


Indices for Spot 6 ■ ( t 2 5 ) 
Indices for Spot 7 » ( t £ 5 ) 





Error rsnge allowed « 3.00 

deg . 

’ * 



Expt.Ang Thert.Ang Optm.ang 

8 3 4 

- - ■■ ' — 

Diffrnce 

Ofit.©iff 






( 

t 

7) 

1-4.36903 

13.70912 

14.96365 

.65992 

1 . 25454 

( 

s 

3) 

18.74132 

21 .02935 

19.47604 

2.28803 

1 .55331 

< 

s 

4> 

27.94303 

27.41825 

28.93628 

.52478 

1 .51803 

{ 

s 

Si 

20.31370 

19.99171 

20.97113 

.32198 

.97942 

i 

s 

6) 

9.69134 

8.59591 

9.98634 

1 .09543 

1 .39043 

1 

z 

7i 

26.78162 

25 . 98589 

27.65150 

.79573 

1 .66561 

( 

3 

4) 

27.67121 

28.42290 

28.72496 

.75169 

.30206 

( 

3 

Si 

29.80225 

32.00727 

30.91108 

2.20502 

1 .09618 

i 

3 

6i 

26.20103 

28.42290 

27.16849 

2.22187 

1 .25441 

f 

3 

T J 

19.53917 

21 .02935 

20.35385 

1 .49018 

,67550 

c 

4 

Si 

IS. 07019 

14.50506 

13.27177 

1 .61487 

1 .23328 

( 

4 

6i 

23.70345 

25.96589 

24.45646 

2.28244 

1 .52941 

( 

4 

7i 

9,97087 

8.59591 

10.27745 

1 .37497 

1 .68154 

t 

5 

&i 

12.61763 

14.50506 

12.98671 

1 .88743 

1 .51834 

( 

5 

7i 

19.86813 

19.99171 

20.5266? 

. 1 2358 

.53496 

( 

6 

7i 

26.62463 

27.41825 

27 . 58488 

.79362 

.16663 


Largest difference between the Angles “ 1.8656 

Smallest difference between the Angles « .0127 

Spot$ 1 S and 

4 belong to zone 

C T 2 0 3 


Spots 1 3 and 

5 belong to zone 

[110 3 


Spots 1 6 and 

7 belong to zone 

C i 1 0 3 


Angle between 1 

Angle between C 

Angle between C 

T 2 6 3 and C 1 1 

T 2 0 1 and Cz 1 

110 3 and Cz 1 

Theoritical 

03 -71 .536263 

03- 36.855064 

0 3 ■108.391327 

Exprimental 

70.000000 

39.500000 

109.000000 

Optimised distance » 2.8626 

Hean square error in angles with 

Optimised Distance 

- 28.823238 


OUTPUT NO. 9 


Sanpl« Nu»b«r: 


6.4 


5. 7.1990 


Number of points indsxetl ■ 6 


Coordinates of the Points 
(X) (Y) 


-a. IS 

1 .05 

-.05 

-1 .35 

.70 

-a . 60 

3.ao 

-a. 80 

-a. 40 

-3.30 

-1 .30 

-4.30 


The index search is between ( S 5 5) and (555) 




Indices 

for 

Spot 

Indices 

for 

Spot 

Indices 

for 

Spot 

Indices 

for 

Spot 

Indices 

.for 

Spot 

Indices 

for 

Spot 


1 - ( a t 1 ) 

a » ( 3 3 4 ) 

3 ■ c 1 a 3 ) 

4 » { 0 1 T ) 

5 » { 4 1 5 J 

6 ■ ( a 1 5 ) 



Error 

range allowed * 3.00 

deg. 

Joints 

■ "t 

Expt . Anu 

a 

Thert . Ang 

3 

Opts. ang 
4 

* 

Biffrnce Ogt.fiiff 

(3-81 <^-3^ 

fei'' ») 

&6 . £3743 
36.01490* 
46.13616 . 

' «T. ‘66 07,3 

83.61907 

34.67347 

43.75464 

87.36869 

':..34..f,:SSS09;!,' 

£5.14883 S.608S4 J .SSSif 

34.66853 = ' 1, . 141.45,:, 

44 . 67,07.1, V , ,.5;; u;9,f,;i^7 



i 

3 

6) 

i 

3 

4 J 

i 

3 

51 

i 

3 

61 

c 

4 

51 

« 

4 

61 

c 

5 

61 


16.59634 

13.99994 

19.54600 

13.67353 

33.35070 

36.53683 

8.30994 


17.53671 
13.95079 
19.03359 
13.74046 
31 .97440 
35.35863 
1 1 .31363 


16.41047 
13.73789 
19.07730 
13.64897 
33.61637 
36 .01917 
8.16358 


.94048 
1 . 04915 
.53341 
. 13307 
1 .37631 
1 . 17801 
3.90368 


1 . 13634 
.78710 
. 05361 
. 09148 
.64187 
.66035 
3.05005 


Largest difference between the Angles ■ 3.1596 
Smallest difference between the Angles > .0536 


Spots 3 £ and 

1 belong to zone 

E 1 5 3 3 


Spots 4 5 and 

3 belong to zone 

E T T T 3 




Theorit ical 

Exprimental 

Angle between C 

1533 and C T T 

T 3 « 52.909093 

55.000000 


Optimised distance « 3.1700 

Hean square error in angles with Optimised Distance ■ 33.340986 


OUTPUT NO. 10 


APPENDIX. IV 


LIST OP VARIABLES 



LIST OF VARIABLE NAMES 


NAME 

TYPE 

REMARK 

1 . ACCEPT 

: Logical 

Decides whether the 
indices assigned to 
any two spots are 
within the error limit 

Z . ALLOW 

: Logical 

Decides whether an 
index selected for a 
particular spot is 
al lowbale . 

5. Dirrci^j) 

: Real, subscripted 

Direrence between 

EAN6 and TANGLE values 

The counters I and J 
represent the spot 
numbers . 

4. DSTAR 

Real 

Optimised value of the 
distance of seperation 

S. DIST 

: Real 

Distance of seperation 

between the sample and 
the X-Ray film. 

6. EANSUiJ) 

: Real, Subscripted 

Mutual angle between 
the Spots obtained 
from their X and Y 
coordinates and the 
initial distance. The 
counters l,J represent 
the Spot numbers. 

7. EANSLEd.J) 

: Real, Subscripted 

Mutual angle between 
the Spots obtained 
from their X and Y 
coordinates and the 
optimised distance. 

The counters represent 
the Spot numbers. 

e. ERROR 

: Real 

Allowable angular deviation 

9. EXTIN 

: Logical 

Decides whether an 
index satisfies the 
extinction conditions. 

10. FINAL 

: Logical 

Limit within which the 
optimised value of 
distance to be found 
out. (S.5 to 3.5 cm! 

It. FINISH 

: Real 

difference between the 
distances of any two 
consecutive iterations 
during optimisation. 

!S. H<I),K(I),L(n 

; Integer, Subscripted 

Indices of the Spots. 

The counter I represents 
the Spot nwmber. 


13. HDIFF 

: Real 


Largest difference 
between EANQ and 

TANSLE values. 

14. lOATE 

: Integer 


Date of execution. 

15. IH{I),IKa),IL(I> 

Integer, 

Subscripted 

Trial indices selected 
for the spots. The counter 

I represents the Spot number 

16 . IP< I , JJ 

: Integer, 

Subscripted 

Numbers of the Spots 
belonging to one zone. 

The counters I and J 
represent the numbers of 
the Spots belonging to 
ine zone. 

17. J( 1 1 

Integer, 

Subscripted 

Counter used to repres 
-ent the position of 
its index in the list. 

The counter represent the 
Spot number. 

18. Ml < 1 ) ,Ma< n ,M3< 1 ) 

5 Integer, 

Subscripted 

indices of Zone axis. 

The counter I represent 
the zone number. 

1». MAX 

Integer 


Maximum limit or the 

Index search 

ffO. niN 

: Integer 


Minimum limit of the 

Index search 

Cl . MONTH 

: Integer 


Month of execution. 

S£. MZN 

: Integer 


The maximum number of 
the zone to which the 
spot belongs. 

S3. N 

: Integer 


Total Number of Spots 
to be indexed. 

C4. NCMMO 

Integer 


Number of times the 
subroutine ERRO is 
cal led. 

C5. NCXf 

: Integer 


Number of times the 
subroutine EXPANG is 
cal led. 

C6. NINO 

: Integer 


Number of iterations 
carried out in the 
subroutine INDEX. 

57. NLAUE 

: Integer 


Number of times the 
subroutine LAUE is 
called. 

E«. NTEST 

: Integer 


Number of times the 
subroutine TEST is 
called. 

6S. NTEI 

: Integer 


Number of times the 
subroutine TRIANG is 
called. 

30. NTEIAL 

: integer 


Number of times the 
subroutine TRIAL is 
called. 



31 . N2 

: Integer 

Total number of zones. 

32. NZH 

: Integer 

Number of times the 
subroutine ZONEM is 
called . 

33, NZONE 

: Integer 

Number of times the 
subroutine ZONECK is 
called. 

34. NZF 

i Integer 

Number of zone points. 

35. N2SB 

: Integer 

Total number of zones 
to which a spot belong 

3&. RANGE 

: Logical 

Decides whether the 
index selected for a 
spot is within the 
range fixed. 

37. RATIO 

: Real 

(C/a) Ratio of Pb-Bi . 

3® . SAMPLE 

: Real 

Sample number. 

39. SDIFF 

: Real 

Smallest difference 
between EAN6 and 

TANGLE values. 

AO. TANGLEtI,J> 

: Real, Subscripted 

Mutual angle between 
the Spots obtained 
from their indices 
using the crystal lo 
-grphic expression 

41 . X< I ) ,Y( 1 1 

: Real, Subscriptd 

X and Y coordinates of 
the Spots on the film. 

42. ZAN6 

: Real 

Mutual angle between 
any two zone axes. 

43. 2AN6LE(I,J) 

: Real, Subscripted 

Mutual angle between 
the zone axes obtained 
after all the spots 
are indexed. 

44. ZERROR 

: Real 

Mean square error 
between EANGLE and 
TANGLE values. 

1 45. ZONANGtlfJI 

: Real, Subscripted 

Mutual angle between 
zone axes, obtained 
from the stereographic 
projection. 

46. ZONE 

: Logical 

Decides whether zone 
checking is required 
for a particular spot. 

47. ZONEC 

: Logical 

Decides whether a set 
of spots belong to the 
same zone. 

48. ZONMAT 

: Logical 

Decides whether the 
mutual angle between 
the zone axes match 


with th# giv*n value*. 




49. 2P(I) 

50. 2SB<I) 


Real« Subscripted 
Integer, Subscripted 


2one point. 

Number of the zone to 
which a spot belongs. 


APPENDIX. 7 


Plow CHARTS 
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